Search templates for stochastic gravitational-wave backgrounds 
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Several earth-based gravitational-wave (GW) detectors are actively pursuing the quest for placing 
observational constraints on models that predict the behavior of a variety of astrophysical and 
cosmological sources. These sources span a wide gamut, ranging from hydrodynamic instabilities in 
neutron stars (such as r-modes) to particle production in the early universe. Signals from a subset 
of these sources are expected to appear in these detectors as stochastic GW backgrounds (SGWBs). 
The detection of these backgrounds will help us in characterizing their sources. Accounting for such 
a background will also be required by some detectors, such as the proposed space-based detector 
LISA, so that they can detect other GW signals. Here, we formulate the problem of constructing a 
bank of search templates that discretely span the parameter space of a generic SGWB. We apply it 
to the specific case of a class of cosmological SGWBs, known as the broken power-law models. We 
derive how the template density varies in their three-dimensional parameter space and show that for 
the LIGO 4km detector pair, with LIGO-I sensitivities, about a few hundred templates will suffice 
to detect such a background while incurring a loss in signal-to-noise ratio of no more than 3%. 
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Multiple earth-based gravitational-wave (GW) detec- 
tors, including both resonant-mass and interferometric 
ones, are currently in operation aiming to make the first 
GW detection. As the sensitivities of these detectors im- 
prove, they will place interesting limits on astrophysi- 
cal event rates and strengths of GW backgrounds, thus 
constraining or falsifying theoretical models. The sub- 
ject of this paper is how to design a template bank for 
searching and bounding the strength of a stochastic GW 
background (SGWB). After formulating the problem for 
a general SGWB, of either astrophysical or cosmologi- 
cal origin, we apply it to the specific case of a SGWB 
with a spectral profile that belongs to a class predicted 
by a host of cosmological models, including inflationary 
and string-theoretic ones. This profile is known as the 
broken-power-law (BPL) spectrum, as described below 

Searching for SGWBs with BPL-type spectra is im- 
portant because some of the cosmological models that 
predict them also allow for their strengths to be large 
enough to be detectable in the near future. In particular, 
in the bandwidth of these detectors, their strengths can 
be several orders of magnitude higher than that predicted 
by the slow-roll inflationary model (SRIM) [I), while 
being consistent with extant observational constraints, 
such as arising from the anisotropy in the cosmic mi- 
crowave background radiation [3|, |^ , the monitoring of 
radio pulses from several stable millisecond pulsars @, 
and the empirical abundances of light elements in the 
universe @. These models are, therefore, especially at- 
tractive to the detectors in operation. Additionally, GWs 
from those unresolved astrophysical sources that have a 
duty cycle appreciably larger than unity will appear as a 



stochastic signal in these detectors, some with BPL-type 
spectra and strengths that may be considerably larger 
than that of a cosmological SGWB predicted by SRIM 
. Indeed in the planned space-based detector LISA |j| , 
the background from unresolved galactic binaries is ex- 
pected to be large enough to form a "source-confusion" 
noise and make it difficult to detect other GW signals • 
Although optimal statistics and templates for searching 
resolvable binaries exist 0, , strategies for character- 
izing such a background in order to allow the detection 
of other sources are beginning to be explored. While we 
derive the limits on template spacings and template num- 
bers for BPL-type cosmological spectra, the formalism 
given here for obtaining the search templates is general 
enough to be applicable to other cases, including astro- 
physical SGWBs. 

We begin by briefly outlining the search statistic for a 
SGWB. The decision on whether a signal is present or 
absent in a detector output is often based on the exam- 
ination of data that are noisy. In decision theory, the 
hypothesis that the data do not contain a signal is called 
the null hypothesis, Hq. Under the alternative hypoth- 
esis, Hi, the detector output is noise plus signal [T^ |. 
Thus, in the frequency domain, the output is 



h(f]0)\ Hl =n(J)+Aw(f;0) 



(1) 



where n(f) is the noise, A is the overall strength of the 
signal and w(f; 1?) gives the spectral characteristics of the 
signal for different choices of the signal parameter vector, 
In general, the waveform will have the appearance: 



i5(/;0) = y(/;0)exp [**(/;<»)] , 



(2) 
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where ?/(/; 1?) is a frequency-dependent amplitude, and 
i?) is the signal phase. We assume that the detector 
noise has a zero-mean Gaussian probability distribution; 
it is described completely by the first two noise moments, 
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n(/) = and n*(/)n(/') = ^P(\f\)S(f - /') , (3) 

where P(\ f\) is the one-sided noise power-spectral density 
(PSD) 13]. We define the inner (or scalar) product of a 
pair of Fourier domain functions a(/) and b(f) as 



(a, bj = 45R f df 



a *(f)kf) 
m(f) 



(4) 



where a(f) and b(f) are the Fourier transforms of tem- 
poral counterparts a(t) and b(t), respectively, and m is 
an inverse weight that, typically, depends on the noise 
PSDs. Its exact form is decided by the detection statis- 
tic at hand. 

The search templates are modeled after the waveform: 



u{f;-d')=Nw{f;V) 



(5) 



where is the template parameter vector and M is a 
normalization factor. The inner-product of a template 
with itself, 



(u(0')» fi(0')) = k(0') 



(6) 



will be taken to be positive definite. Above, *Jk is the 
template norm, which is parameter dependent, in gen- 
eral. The normalization factor is related to the template 
norm as follows: 



-1/2 



(7) 



To test a hypothesis, one computes the cross-correlation 

of the data with the templates, viz., m(i?')Y 

which is termed as the matched-filter output (MFO). Un- 
der Hi, the mean of the MFO is 



(8) 



where S'd = t? — 1? . 

One often uses in searches unit-norm templates, 
namely, 



k(0') 



1 . 



(9) 



The advantage of using such templates is that under Hi 
and for S'd = 0, the mean of the matched-filter output 
(MFO) is just the signal strength divided by the template 
normalization factor, i.e., 



S{djd = 0) := [h(0)\ Bl , «(*) 



A 



7V(#) 

assuming that the signal model is perfect. 



(10) 



To quantify the effect of a mismatch, S'd, it is useful 
to introduce the match or ambiguity function: 



M(0,<W) := (u(0), fi(0')) , 



(11) 



which tends to k(#) as S'd — > 0. Then the mean of the 
MFO can be shown from Eq. © to be 
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M (0, fa?) . 



(12) 



For small values of S'd, one can Taylor expand M about 
S'd = to obtain 

M(0,6#) = K(-d) \l + 

-g^WS&'SV] + 0(6# 3 ) (13) 

where the Einstein summation convention over repeated 
indices, \i and v, was used and we defined 



da := - 

K 



dM 



d8d» 



S-d=0 



1 

2~- 



d 2 M 



dSd^dSd" 



<W=0 

Above, t/ M „ can be interpreted as the metric on the pa- 
rameter space that maps parameter mismatches into dips 
in the signal-to-noise ratio (SNR) [lj| , provided van- 
ishes. (The MFO of a unit-norm template is equivalent 
to the SNR 0.) 

It is important to note here that an observer also has 
the choice of using unnormalized templates, such that 
J\f = 1 in Eq. (JjJJ. This has the advantage that one does 
not have to recompute M and, therefore, the search tem- 
plates, for every value of i? or every time m (which can 
be the noise PSDs of the detectors) changes. Indeed, this 
choice was exercised by Ungarelli and Vecchio in their pi- 
oneering work in Ref. Q. However, the disadvantage of 
such a choice is that the associated ambiguity function, 
M(i?, S'd), has first order errors arising from parameter 
mismatches. Consequently, a "wrong" template (i.e., a 
template with S'd ^ 0) applied to a given data set can 
actually trigger an MFO that is larger than that of the 
"correct" template (with 6$ = 0) applied on the same 
data set. Use of constant-norm templates, such as the 
unit-norm ones defined above, avoids that problem. 

To see explicitly why d M need not be zero for unnor- 
malized templates, note that the associated ambiguity 
function obeys the Cauchy-Schwarz inequality [16(: 



M{&,8&) = 2(0')) ^ J *(#)*(■&') . (I 5 ) 

Since, in general, ) can be larger than for some 
7^ -d, the rhs above can actually exceed k(i?). Thus, 
need not be the maximum value of M(t9, S'd) for un- 
normalized templates. Equation 1)13(1 then implies that 
d^ need not be zero for such templates. However, for 
unit-norm templates the rhs of Eq. I|15fl is identically 
unity (and k(i9) = 1 = independent of the value 

of i9' or i9. There, M("d, S'd) attains the maximum pos- 
sible value of 1, when 'd' — -d. Thus, d M has to vanish for 
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unit-norm (and constant-norm) templates, and can 
assume its role as a parameter-space map. 

Without any means for distinguishing a stochastic 
GW background in a detector from the detector's in- 
trinsic noise, the search for such a signal involves cross- 
correlating the outputs of a pair of detectors. As shown 
in Ref. |l7j . a useful statistic in decision making in this 
context is the cross-correlation (CC) statistic, 



S 



df 



df'5 T (f'-f)h A (f')h B (f)u(f) , (16) 



where h A (f) is the inverse Fourier transform of the strain 
in the Ath detector, T is the observation time, Srif) = 
/— t/2 ^ exp(— i27r/t) is the finite-time approximation of 
the Dirac delta function, and u(f) is a filtering function 
that will be determined below. The CC statistic can also 
be cast as the output of a matched filter: 



S= K 



\h x , u 



(17) 



where K\h x is a functional of a pair of detector inputs: 



K(f) := - 



df 



,h A {f)h* B {f)W - I) 
PaU)PbU) 



(18) 



and the inner product is defined as in Eq. with the 

[P A (f)P B (f)]-\ 
appearing in 5* is a ran- 



inverse weight there set to m = 

The product \h* A {f) x h B (f) 

dom variable since the SGWB strains, produced cos- 
mologically or astrophysically (in some cases), are so. 
The detection statistic, therefore, is the mean of the CC 
statistic. 



s={k\ Hi , 

And the variance of 5* is 



(19) 



(20) 



which defines the noise-squared of S. If one assumes that 
the noise power in each detector due to terrestrial sources 
is much larger than that due to a SGWB, then 01 



which is maximized when the filtering function matches 
the signal, i.e., 



u = \K 



(23) 



where A is a proportionality constant. Although S and 
a are dependent on the choice of this constant, the SNR 
itself is independent of it. 

Equation \'1'M suggests u as templates for searching 
an astrophysical or cosmological SGWB, as long as the 
assumptions made above remain valid. We now concen- 
trate in the rest of the paper on the search templates re- 
quired for a cosmological SGWB. (A similar problem for 
the astrophysical SGWB will be studied elsewhere [IS}-) 
Theoretically, the strain due to a cosmological SGWB 
in each detector is expected to have a Gaussian probabil- 
ity distribution with zero mean; their variance-covariance 
matrix elements are given by [ltj: 



h* A (f) 



M/0 

Hi Hi 



20tt 2 

xS(f-f') 



|/|- 3 7ab(|/|)0gw(|/|;^) 

(24) 



where Hq is the Hubble constant, ^y AB is the overlap- 
reduction function (ORF) for the detector pair [2(]], 0gw 
is the energy density of the stochastic GW per logarith- 
mic frequency bin divided by the critical energy density 
required to close the universe, and # are the signal pa- 
rameters on which it depends. Heretofore, we will iden- 
tify j AB = 7. The ORF for co-located and co-aligned 
interferometric detectors with orthogonal arms is nor- 
malized so that it is identically unity. Using the above 

strain-power density in the expression for K\ Bl yields the 
cosmological SGWB template: 



= N 



7(I/I)0gw(I/I;0) 



\f\ 3 P A (\f\)P B (\f\) 

where, for unit-norm templates, 



N = 



df 



7 2 (l/l)0 2 G w(l/l;^) 

pp A (\f\)p B (\f\) 



-1/2 



(25) 



(26) 



We now show how to obtain the spacing between such 
templates on the 1? parameter space such that its dis- 
creteness, Sfl, is small enough to guarantee an SNR of 
>97% of that obtained in the ideal case of S& = 0, 



— T 

T 



df \u(f)\ 2 P A (\f\)P B (\f\) 



(21) 



Above, it was assumed in the second approximation that 
the cross-correlation of the terrestrial noises in the two 
detectors is negligible. Thus, the SNR is 



(22) 




1. Single power-laws 

Before studying the template spacings of the BPL spec- 
tra, let us consider the simpler case where OqwC/) is m 
the form of a single power-law (SPL) in frequency, 



Ogw(/) =^0(7- 
Jp 



(27) 



where k is a real power, and Sl an( i fp are positive- 
definite real constants. In the bandwidths of LISA or 
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LIGO-type detectors, the SGWB spectrum predicted by 
SRIM will appear as a special case of the above, with k w 
0. In such a case the only intrinsic search parameter is 
the index k and the ambiguity function can be expanded 
around 5k := k — k' = as 

M(k,5k) = l-^[a(2)-a 2 (l)]5k 2 +0(5k 3 ) , (28) 

where a (n) is the template-space mean of In (f/f p ): 

a(n) = (fi(fc), [In (f/f p )] n u(k)) 

f + /f » dx l 2 {f p x)x 2k - Q [\nx] n 



j f. P A {fpX)P B (fpX) 
f+/f, ^(y^gM-e 
/-//, P A{fpX)P B {fpX) 



-i -1 



(29) 



is readily deduced from Eq. I|28|l as 



9ll (k) = -[a(2)-a 2 (l)] 



(30) 



which has the following properties: First, the Cauchy- 
Schwarz inequality can be used to prove that gn(k) is 
non-negative for all real k Second, Eqs. Qg) and 

(|29|l show that it is dependent on k, which is confirmed 
by Fig. ^ This implies that for the optimal coverage of 
the A;-space the template-spacings must be chosen to vary 
in step with the gn(k) values. The template spacing is a 
minimum at the global maximum of gu(k) (at k ~ 0.5), 
and is: 



elk = 0.637 



1-MM 
0.01 



1/2 



(31) 



Above, n is any real number, and /_, /+ are the lower 
and upper limits of the frequency integral, respectively. 
For the LIGO-Hanford (LHO) and -Livingston (LLO) 
pair, with LIGO-I sensitivity, we choose /_ = 40Hz 
(determined by the detectors' seismic-noise cut-off) and 
/+ = 512Hz, respectively. For this inter-site correlation, 
the statistic does not receive any appreciable contribu- 
tion from higher frequencies owing to a weak "f{f). By 
contrast, for co-located interferometers, j(f) — 1 for all 
frequencies, and the upper cut-off frequency is deter- 
mined by the worsening sensitivity due to the photon 
shot-noise |0, |2^. An optimal way of computing the 
CC statistic on the data from a set of multiple detectors 
that includes a co-located pair, while accounting for 
intra-site terrestrial noise, was obtained in Ref. [2l| . 
For searches in this kind of a detector network, the 
value of f + should be revised upward of 512Hz. As 
such, the expressions here are applicable to any pair of 
GW detectors. However, the template-spacings, number 
of templates, and the figures are computed for the 
LLO-LHO (4km) pair, with LIGO-I sensitivities El . 




where MM is the minimal match required by an observer 
between the discretely spaced templates and a signal. 
Typically, MM is set equal to 97%. But we find above 
that dk is as large as 0.637 for a minimal match as high 
as 99%, as is evident in Fig. [21 




Power-law exponent k 



FIG. 2: The ambiguity function, M(k, 5k), for a single power- 
law SGWB plotted as a function of k and the mismatch, 5k. 
For any given k, the function attains the maximum possible 
value of unity when 5k — 0. And for any given 5k, the function 
is a minimum at k = 0.5000tQ 0625 ' which is consistent with 
the behavior of the metric depicted in Fig. HI |2ij|. 



If one were to choose the template-spacing to be uni- 
formly equal to the above value and, therefore, err on the 
side of over-covering the parameter space, then the num- 
ber of templates required for a search with k € [—4, 4] 
and a minimal match of at least 99% is I24JI 



FIG. 1: The only metric component, gn, for a single power- 
law SGWB plotted as a function of k. The template spacing 
is a minimum near the global maximum at k = 0.50001 



-,+0.0625 
-0.0625 • 



The only metric component for the above SGWB signal 



N = Ceiling 



Mnax "'mm 



0.637 



13 



(32) 



which is easily implementable in real time on the data of 
a detector pair. 
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2. Broken power-laws 

A likely character of 51gw(/) will be in the form of a 
broken power-law (BPL) in frequency |2|, 



Ogw(/) = f2o 



e l- 



+ e 



J_ 

f P 



- i 



(33) 

where fc_ > 0, fc+ < are real power-law exponents, 
and the peak frequency, f p , and f^o are positive-definite 
real constants. The first three are intrinsic search pa- 
rameters, which define the three-dimensional parameter 
space, "d = (f p /fpo, k-, Here f p o is an arbitrary 

reference frequency chosen large enough so that 5f p /f p o 
is small and the O (<5# 3 ) terms in Eq. (|13[) are indeed 
negligible. 

We calculate the ambiguity function for the templates 
in Eq. I|25|) . with the above ^gw(/)) and compute the 
metric from its second derivatives using Eq. I|14|) . These 
derivatives now involve the following integrals of the log- 
arithm and different powers of /: 



D 



dx: 



l 2 (f P x)x 



2fe_-6 



f-/fp 



PaU p x)PbU p x) 



f+/fp 



-I -1 



dx 



a± (n) 

0±(n) 
(3(n) 



Y(f p x)x 2k + 
P A {f P x)P B {f p x) _ 

D l dx ^(f P x)x^- e [Ins] 

>f±/f P PA(f p x)P B (fpX) 

I" 1 f ^\ry,2k±~6 



:= TP>k 



dx 



rU P x)x 2 



f ± /f p PA{fpX)P B (fpX) 

/3_(n)+/3+(n) , (34) 



where the dependence of D, a±(n), (3±(n), and (3(n) on 
the detector indices A and B is implicit. 

The six independent components of the symmet- 
ric metric on the three-dimensional parameter space, 
(f P / fpO,k-,k+), are: 





0.4p 




FIG. 3: The template-density p (normalized here to have a 
maximum value of 1) for a BPL parameter-space plotted as 
a function of k± and f p , two variables at a time. The fixed 
parameters in these plots are (clockwise from top) f p = 50Hz, 
k- = 0.125 and k+ = -0.125. 




9oo 
9n 

922 

9oi 

902 
912 



i(^) [/3(2)-/3 2 (l) 
\[a + (2)-al(l)) , 



fpQ 
U 

fpQ 
f P 



fc_]a_(l) 
k+]a+(l) 



1 



--a_(l)a+(l) , 



(35) 



where the parameter vector components were taken to 
be $° = fp/fpo, A 1 = k~, and d 2 = k + . Once again it 



FIG. 4: The ambiguity function, M(i9,<5#), for a broken 
power-law SGWB plotted as a function of f p and the mis- 
match, Sf p (both in Hz). Here, we set k± = =p2 and 5k± = 0. 
Note how quickly M asymptotes to a value of unity for large 
f p , implying a small number of templates for f p ~ 200Hz. The 
ambiguity function is within 99% for: (a) f p £ [50, 120]Hz, if 
Sfp < 4Hz, (b) f p 6 (120, 150]Hz, if Sf p < 10Hz, and (c) 
f p > 150Hz, if 6f p < 20Hz. 



is clear that the metric is dependent on the parameters, 
as confirmed by Fig. [3J Thus, the template spacing for 
a fixed minimal match is not uniform |25j |. 

The number of templates, N, can be estimated from 
the above metric by dividing the parameter space volume, 



FIG. 5: M(6, 86) for BPL spectra plotted as a function of k+ 
and k- for four values of f v = 50, 100, 190, 450Hz (clockwise 
from top left). The mismatch values are Sk± = 0.25 and 
Sf p — 2.5Hz. As f p increases, note how little M varies with k+ 
(for any given fc_ ) . This is indicative of the template spacings 
being larger along k+ than along fc_. This is expected since 
most of the contribution to the CC statistic arises at lower 
frequencies, where the fc_ index has support (especially, when 

f p ~ 200Hz). 



V = / /det 



(36) 



by the volume of a unit cell in three-dimensional space, 
v = (2yJ(l - MM)/3 ) 3 0. We numerically compute 
the above volume and find that 



N ~ 250 



1 



0.03 



MM\ 



-3/2 



(37) 



FIG. 6: M{6, 86) for a broken power-law SGWB plotted as 
a function of k- and f p for four values of fc+ = -4, -2, -1, 

The mismatch values are 8k± = 



(clockwise from top left). 
0.25 and 8f p = 2.5Hz. 



. 9951 
M 0.991 
. 985, 
. 98" 
0. 97" 



where we again err on the side of over-coverage by re- 
laxing the numerical computation to allow for MM > 
97%. The parameter-dependence of the template den- 
sity, p — ^/det Wg^vW/v, is illustrated in Fig. 01 The 
above value of N is consistent with the near-unity value of 
the ambiguity function shown in Figs. 14171 for template- 
spacings as large as 5k± = 0.25 and Sf p = 2.5Hz. 

It has been projected in Ref. [l7J that LIGO-I and Ad- 
vanced LIGO may succeed in placing upper limits on CIq 
of the order of 5 x 10 -6 and 5 x 10~ n , respectively, for 
the k = SPL spectrum. The first science run at LIGO 
already demonstrated successfully the application of a 
single template (i.e., the k — case of SPL) on the data 
from the LIGO detector pairs to obtain bounds on £1 
[2^,123. With the upcoming science runs at LIGO, the 
sensitivities are fast approaching closer to the designed 
target so as to make the first upper limit given above 






450 M 




FIG. 7: M(6, 86) for a broken power-law SGWB plotted as 
a function of k+ and f p for four values of k- = 0, 1, 2, 4 
(clockwise from top left). The mismatch values are 8k± = 
0.25 and 8f p = 2.5Hz. 



realizable in the near future. This progress necessitates 
the availability of techniques and templates to look for a 
variety of proposed astrophysical and cosmological SG- 
WBs in the ever-so sensitive data. This paper addresses 
this issue for the latter category of signals, which assumes 
the background to be isotropic and unpolarized. The for- 
mer case of an astrophysical background will be discussed 
elsewhere [l8| . 
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